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As  sound  pressure  increeises  to  "finite"  levels » nonlinear  effects  usually 
ignored  in  acoustics  problems  become  increasingly  Important.  In  this 
paper  plane  waves  propagating  through  a pipe  are  considered.  In  the 
first  section  the  distortion  of  a wave  that  is  sinusoidal  at  the  source 
is  treated.  In  the  second  section  the  source  signal  is  a pulse  of  band- 
limited  (500  to  3500  Hz)  random  noise.  In  both  cases  theoretical  results 
are  compared  with  experimental  measurements.  Extension  of  the  analysis 
to  the  more  practical,  case  of  outdoor  propagation  is  indicated. 
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A good  deal  of  fundamental  research*  both  experimental  and 
theoretical*  has  been  done  on  finite  amplitude*  waves.  Almost  all  of  it 
has  been  done  on  deterministic  waveforms;  very  little  has  been  done  on 

(l  2) 

finite  amplitude  noise.  * The  ultimate  goad,  of  the  present  research 
is  to  attack  the  problem  of  propagation  of  finite-amplitude  noise.  In  an 
effort  to  simplify  ths  problem  as  much  as  possible*  it  was  decided  to 
restrict  the  investigation  to  pleme  waves.  But  how  does  one  handle  a 
finite  amplitude  wave  whose  waveform  is  not  deterministic?  Our  approach 
has  been  to  use  weak -shock  theory*  cast  in  the  form  of  a computer 
algorithm,  to  obtain  theoretical  predictions  and  to  make  measurements  in 
a plane  wave  tube  to  check  these  predictions.  To  test  this  approach  we 
have  made  a preliminary  study  of  initially  sinusoidal  waves.  Thic  has 
been  followed  by  a study  of  the  propagation  of  pulses  of  random  noise. 

The  noise  is  Gaussian  and  bandlimited  to  approximately  3 octaves.  Sound 
pressure  levels  up  to  l60  dB  (re  0.0002  wbar)  have  been  used.  The 
results  of  these  studies  form  the  basis  of  this  paper. 

Use  of  the  plane  wave  tube  has  introduced  the  complication  of  tube 
W8dl  attenuation  and  dispersion.  Accordingly,  ve  have  had  to  modify  weedt- 
shock  theory  to  take  account  of  these  effects.  It  is  simpler  in  the 
beginning,  however,  to  consider  distortion  in  the  absence  of  tube  wall 
effects . 


"Finite  amplitude"  is  a qualitative  terra.  In  genereJ.  it  implies  that 

nonlinear  effects,  which  are  ignored  in  linear  acoustics,  must  be 

accounted  for.  However,  some  distance  may  be  required  for  the  nonlinear 

effects  to  be  noticeable.  For  example*  in  air  an  initially  sinusoidal 

plane  wave  of  SPL  130  dB  (re  0.0002  ybar)  forms  shocks  at  a distance  of 

about  210  wavelengths.  An  observer  much  closer  than  this  to  the  source 

would  probably  view  the  wave  as  "infinitesimal"  whereas  an  observer  much 

further  away  would  see  clear  evidence  that  the  wave  "finite  , 

amplitude."  V 3 
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When  sound  waves  of  large  aaplitude  propagate,  the  nonlinearity  of 
the  prfjpaeHtion  process  causes  progressive  distortion  of  the  waveform. 

In  a nondispersive  medium  this  distortion  will  cause  a sine  wave  to 
deform  into  a sawtooth  wave,  provided,  of  course,  the  wave  does  not  first 
become  ^’infinitesimal"  through  absorption.  Figure  1 depicts  this  pro- 
gressive distortion  for  an  arbitrarily  shaped  wave.  The  vertical  axis  is 
particle  velocity  (or  pressure)  and  the  time  base  is  retarded  time 
t’*t-x/c^.  Note  how  the  portions  of  the  waveform  having  larger  paurtlcle 

velocities  tend  to  overtake  the  slower  velocity  poi^icxis.  When  this 
overtaking  actually  occurs,  shocks  form.  This  is  clearly  shown  at  x«x. 
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FIGURE  1 

WEAK -SHOCK  THEORY 


Mathematically,  the  progressive  distortion  of  the  continuous  portion 
of  the  waveform  ir  described  by  the  propagation  lav; 


Here  x is  distance,  t time,  u particle  velocity  and  c small-signal  sound 

speed;  for  gases  6 = (y  is  in  the  ratio  of  specific  heats)  and  for 

liqiiids  3 = ^ (b/A  is  a nonlinearity  parameter).  At  the  discontinuity 

the  propagation  of  the  shock  front  is  described  by  the  weak-shock 
relation 
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where  and  are  values  oi'  u Just  ahead  of  and  Just  behind  the  shock, 

respectively.  Notice  how  different  points  on  the  waveform  move  with 
respect  to  each  other.  As  dissipation  at  the  shock  front  begins  to  take 
its  toll  of  the  wave,  points  of  the  waveform  in  the  neighborhood  of  the 
shock  Just  disappear.  In  the  end  only  those  sections  of  the  waveform 
that  were  originally  near  the  zero  crossings  survive. 

Althou^  an  explicit  analytical  solution  of  Eqs.  (1)  and  (2)  (and 
appropriate  auxiliary  equation)  can  be  found  for  simple  waveforms,  such  a 
solution  is  not  feasible  for  c<»aplicated  waveforms,  particularly  non- 
deterministic  ones  such  ais  noise.  A computer  ailgorithiB  in  which  Eqs.  (l) 
and  (2)  are  used  as  the  basis  for  calculating  the  distortion  of  an 
arbitrary  waveform  is,  however,  possible. 

We  have  developed  such  an  algorithm.  Equations  (l)  and  (2)  are  cast 
in  the  form  of  difference  equations  and  solved.  I^e  wave  is  "cos^uter 
propagated*'  a small  distance  x,  and  the  original  time  base  is  distorted 
in  accordance  with  the  relation 


t' 

new 


^old 
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bid 


x/c. 


(3) 


This  distorted  time  base  is  scanned  for  multivaluedness.  If  none  is 

found,  another  small  distance  step  is  taken.  'Ibis  sequence  is  repeated 

until  finallj*  multivaluedness,  which  signifies  the  presence  of  a shock,  is 

found.  The  shock  is  located  as  prescribed  by  Eq.  (2)  and  the  particle 

velocities  at  the  shock  front  are  adjusted  accordingly.  In  a periodic 

waveform,  such  as  a sawtooth  wave,  the  shock  fronts  are  stationary  in  the 

i*etorded  time  frame  because  the  shocks  move  with  speed  c . However,  in 

o 

an  arbitrary  wave  field  the  shocks  exhibit  relative  motion.  A separate 
portion  of  the  distortion  program  monitors  this  motion  and  accounts  for 
shock  merging.  Ibrough  the  progression  of  growing,  decaying,  and  merging, 
the  shocks  ultimately  determine  the  shape  and  amplitude  of  the  wave.  The 
distortion  algorithm  has  been  tried  on  several  initial  signals,  including 
N waves , sinusoids , and  rcuidomly  generated  waves . The  results  are  in 
excellent  agreement  with  such  analytical  results  as  are  known. 

Unfortunately,  this  simple  algorithm  does  not  give  a very  good 
account  of  itself  when  the  wave  motion  is  in  a tube.  In  19615 » McKittriek, 

(3) 

Blackstock,  end  Wright  reported  on  the  measurement  of  the  wave  shape 
of  repeated  shocks  in  a tube.  Rather  than  exhibiting  the  sawtooth  shape 
that  was  expected,  these  waves  had  a ma^'kedly  rounded  top  while  retaining 
a sharp  trough.  It  was  demonstrated  thi  t the  asymmetry  was  due  to  dis- 

(U) 

persion  caused  by  tube  wall  effects.  In  a theoretical  study  Coppens 
confirmed  the  asymmetry  and  its  cause  (dispersion)  for  waves  in  the  pre- 
shock region.  We  have  also  observed  the  same  asymmetric  waveform  in  our 
measurements.  It  should  be  noted  that  the  tube  wall  caiises  attenuation 
Eis  well  as  dispersion. 


Our  Method  of  accountiaf;  for  tube  uull  effects  is  to  incorporate 
attenuation  and  dispersion  directly  in  our  distortion  algorithm.*^  See 
Fig.  2.  Briefly,  for  each  incresental  propagation  step  the  vave  is 
^ppropriately  distorted.  *n)en  an  FFT  routine  is  used  to  express  the  vave 
in  terra  of  its  spectrsl  coaponents.  The  coaplex  as^litude  of  each 

spectral  component  is  corrected  for  the  attenuation  and  phase  shift  that 
occurs  over  the  increnental  step  as  follows: 


C ' 
n 


ox(l+J) 


(U) 


where  C ' is  the  corrected  aaplitude  and  a is  the  attenuation  coefficient 
n 

for  the  fundaaental.  As  is  well  known  for  tube  wall  attenuation  a is 

proportional  to  /w,  and  this  is  ^rtiy  the  factor  /n  appears  for  the  n^" 
harsionic.  Note  also  the  phase  shift,  which  is  caused  by  the  dispersion. 
Since  the  phase  shift  is  different  for  each  harnonic,  the  spectral  cob> 
ponents  get  sli^tly  out  of  step  with  each  other,  this  is  what  causes 
the  syBBsetry  of  the  vavefora  to  be  ruined,  particularly  at  the  shocks. 

The  asyaratry  has  an  effect  cxi  the  overall  decay  of  the  vave,  as  will 
presently  be  explained.  The  waveform  is  inverse  transformed  and  another 
distance  step  is  taken.  In  this  way,  the  wave  profile  at  any  distance 
from  the  source  is  realised. 
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FIGURE  2 

A SCHEMATIC  DIAGRAM  OF 
MOOI'^IED  WEAK-SHOCK  THEORY 

Now  let  us  txirn  to  the  experiment.  The  plane  wave  tube  is  a 96  ft, 

2 in.  i.d.  aluminum  pipe  cosposed  of  eight  12  ft  pipe  sections.  The 
sections  are  Joined  by  carefully  constructed  2-piece  flanges  that  are 
designed  to  minimize  reflections  at  the  Joints.  The  pipe  is  tenrlnated 
in  a suitable  nonref looting  fiberglass  wedge.  Source  frequencies  in  the 
range  800  to  3500  kHz  and  SPL's  in  the  range  lUo  to  l60  dB  have  been  used. 
In  the  following  results  the  fundamental  frequency  was  2 kHz  and  the 
initial  SPL  was  l60  dB.  The  sound  source  was  an  ALTEC  290  E "Giant  Voice" 

• ' («;) 

The  algorithm  thus  becomes  similar  to  one  used  by  Cook  to  calculate 

the  distortion  of  a sinusoidal  vave  in  a free  thermoviscoxis  medium. 

Cook  did  not,  however,  account  for  shocks  in  the  waveform  by  means  of 

weak-shock  theory. 


horn  driver  powered  in  a pulsed  node  by  s 200  W Duksne  power  SKplifier. 
The  pickup  elcrophone  vss  a B k K 1/U  in.  type  U136,  positioned  in  a 
aount  so  that  the  diaphra«n  was  flush  with  the  pipe  wall . The  microphwie 
output  was  displayed  on  an  oscilloscope  and  the  wavefora  i^otographed. 
Figure  3 shows  a generalized  diagrait  of  the  equipment. 
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FIGUftE  3 

A GENERALIZED  EQUIPMENT  ARRANGEMENT 


Figure  U shows  one  cycle  of  the  measured  and  ccaaputed  wavefonas  at 
distences  of  1 ft  and  13  ft  from  the  sound  source.  Ihe  waveform  at  one 
foot  is  already  scaaewhat  distorted  primarily  because  of  nonlinearities  in 
our  sound  so’orce.  This  vaveform  is  the  input  signal  for  both  the  experi- 
ment and  the  computations.  At  13  ft  the  wave  has  traveled  about  3 'hock 
formatitvi  lengths  so  that  on  the  basis  of  simple  weak-shock  theory  it 
would  be  expected  to  be  essentially  perfect  sawtooth.  Because  of  wall 
effects,  howe'/er,  there  is  already  some  rounding  at  its  peak.  Notice 
that  only  the  peak  Is  rounded,  not  both  the  peak  and  the  trough.  As 
explained  earlier,  the  asymmetry  is  due  to  dispersian.  The  dashed  wave- 
form was  calculated  by  omitting  the  wall  effects  portion  of  the  algorithm. 
Thus,  this  is  the  waveform  one  would  expect  in  an  open  medium  and,  of 
course,  is  the  same  one  could  cwopute  analytically  from  ordinary  weak- 
shock  theory. 
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In  Tin,  5 13  ^ co^tikrlson  been  r«p«nt«d  and  •xp«riiMntaI 

•nd  cot^uted  vnwforas  nt  ft  and  85  t%  added.  Mote  the  tncreaeed 
aayMetry»  particularly  the  roundinti  of  the  peak , as  the  wave  propagatet. 
Although  the  wave  depicted  here  ia  a relatively  strong  <»e,  we  have  used 
our  Mthod  to  treat  weaker  waves  with  equally  good  results.  The  dif- 
ference between  the  s olid  and  dashed  line  theoretical  curves  shows  how 
ii^rtar.t  tube  wall  effects  arc  under  our  experiaental  conditions. 

In  Pig.  6 three  theoretical  curves  of  the  decoy  of  the  wave  with 
distance  are  shown  together  with  our  cKperiaental  results.  The  top  curve 
is  for  ordinary  weak-shock  theory.  Again » this  is  the  sane  at  our  co«- 
puter  Method  gives  if  tlte  effects  of  the  tube  wall  are  onitted  entirely. 
The  only  attenuation  la  that  associated  with  the  dissipation  taking  place 
at  the  shocks.  The  middle  curve  is  obtained  froo  our  cooqputcd  waveforas 
and  thus  represents  a model  in  which  t*.be  wall  attenuation  and  dispersion 
are  included  along  with  shock  dissipation.  The  lowest  curve  is  for  weak- 
shock  theory  corrected  for  tube  wall  attenuation  but  not  for  dispersion. 
To  obtain  this  cxirve  we  have  uced  the  equation  first  suggested  by 

Rudnick,^^^ 


’ip, 


max 


dx 


- up 


max 


wpo  W 


(5) 


This  equation  gives  the  rate  of  decay  of  a sawtooth  wave  as  the  sum  of 
the  deesy  due  to  ordinary  attenuation  (the  first  tern)  and  the  decay  due  to 
dlssipatlwi  at  the  shocks  (the  second  tern).  Here  p^^^  is  pesk  pressure, 

X distance,  c«u  /c  , k wave  nunber,  p the  peak  pressure  at  x«0,  snd  & as 
0 c o 

defined  previously.  The  value  of  a is  given  by  Rudnick's  summation  method 
for  calculating  the  effective  tube  wall  attenuation  coefficient. 


The  experine.ital  points  confirm  the  middle  curve.  The  top  curve 
represents  a pi^dictiwj  of  too  little  attenuatior.  because  tube  vail 
absorption  Is  neglected.  Y«t  when  wall  absorption  is  accounted  for,  one 
obtains  the  lowest  curve,  and  too  nuch  attenuation  is  predicted.  We 
believe  that  the  key  to  the  success  of  our  model,  i.e.,  the  middle  cu.<e, 
the  dispersion.  The  severe  rounding  of  the  peaks  cf  the  shocks  caused 
by  dispersion  in  effect  reduces  the  amplitude  of  each  shock.  Instead  of 
extending  from  trou^  to  peak,  the  actual  shock  runs  only  froo  trough  to 
a point  short  of  the  peak.  This  lessens  the  shock  dissipation.  In  other 
words,  dispersion  inhibits  the  attenuation  associated  with  the  shock. 

When  weak-shock  theory  is  modified  to  take  account  of  both  dispersion 
and  attenuation,  the  result  is  the  middle  curve. 


The  decay  of  sawtooth  waves  in  tubes  has  received  the  attention  of 

several  previous  investigators . In  general  they  found  a discrep- 

ancy between  their  measured  values  of  decay  auid  predicted  values  obtained 
by  adding  wall  losses  to  shock  dissipation.  As  we  now  know,  dispersion 
should  have  been  taken  into  accovint  as  well  as  wall  attenuation.  However, 
when  we  compared  their  data  with  predictions  based  on  our  theory,  in 
which  account  is  taken  of  dispersion,  completely  satisfactory  agreement 


sy/ 
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FIGURE  6 
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was  still  lacking.  The  reasons  for  this  are  not  clear,  in  that  we  have 
found  good  agreement  with  our  owTi  measurements  for  a variety  of  frequen- 
cies in  the  range  830  to  3500  Hz  and  sound  pressure  levels  in  the  range 
lUO  to  160  dB.  There  ai*e  some  differences  between  the  earlier  measure- 
ments and  o\irs.  We  have  generally  used  higher  frequencies  and  somewhat 
Icwer  sound  pressure  levels.  Further,  the  emphasis  in  the  earlier  work 
was  on  peak-to-peak  amplitude.  Little  attention  was  paid  to  wave  shape 
so  that  we  do  not  know  the  extent  to  which  rounding  of  the  peaks  of  the 
waveform  (Figs.  U and  5)  might  have  influenced  interpretation  of  the  data. 
Under  the  circijmstances , all  that  can  be  said  definitely  at  this  point  is 
that  our  theoretical  method  has  been  verified  for  the  middle  audio  ranges 
for  sound  pressure  levels  up  to  about  I60  dB. 


In  order  to  use  the  previously  described  algorithm  for  noise  pulses , 
only  a few  program  changes  are  necessary.  The  fundamemtal  frequency 
becomes 
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fund 


1/T 


(6) 


where  T is  the  length  of  the  pulse.  Attenuation  and  dispersion  are  again 
applied  to  each  spectral  component  up  to  the  Nyquist  frequency. 

The  experimental  measurements  were  made  on  a perticuletr  pulse  of 
noise  as  follows:  First  a random  electrical  pulse  of  noise  was  bandpass 

filtered  and  recorded  on  a uape  loop.  The  pulse  was  then  played  back, 
amplified,  and  applied  to  the  sound  source,  in  this  case  a University 
65  W horn  driver.  As  a result  of  the  filtering  action  of  the  electrical 
filter  and  the  horn  driver,  the  spectrum  of  the  acoustic  noise  pulse  as 
it  started  down  the  tube  extended  from  about  500  Hz  to  about  3500  Hz. 

The  initial  (pulse)  SPL  was  I60  dB,  and  the  pulse  repetition  rate  was  one 
per  second.  A memory  type  oscilloscope  was  used  to  record  the  waveform 
of  the  p\ilse  at  various  distances  down  the  pipe. 

In  Fig.  T the  experimental  and  computed  waveforms  are  compared. 

Note  how  the  noise  signal  simplifies  as  the  wave  propagates. 

The  number  of  zero  crossings  diminishes  and  minor  irregularities  in 
the  waveform  disappear  as  they  are  "eaten  up"  by  the  larger  shocks.  Note 
also  the  roimding  of  the  wave  peaks.  As  in  the  sinusoidal  case,  this  is 
caused  by  dispersion.  It  will  be  seen  that  the  agr^^ment  between  theory 
and  experiment  is  quite  good. 

In  terms  of  frc'^uency,  the  reduction  in  number  of  zero  crossings 

represents  a transfer  of  energy  to  xhe  lower  end  of  the  spectrum.  On  the 

other  hand,  the  presence  of  the  shocks  with  their  very  fast  rise  times 
means  that  the  high  frequency  end  of  the  spectriim  has  also  been  enhanced. 
In  general  the  relative  buildup  of  the  upper  and  lower  ends  of  the  spec- 
trum takes  place  at  the  expense  of  the  middle  of  the  spectrum.  The 

effect  of  nonlinearity  is  thus  to  flatten  out  the  spectrum. 

In  Fig.  8,  the  13  ft  and  05  ft  computed  waveforms  from  Fig.  T are 
repeated  and  comprored  to  waveforms  computed  by  using  the  algorithm  based 
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on  weak-shock  theory  alone.  In  the  latter  case,  no  correction  was  made 
for  attenuation  and  dispersion  caused  by  wall  effects.  It  is  clear  that 
the  wall  effects  have  a significant  influence  on  the  waveform,  particu- 
Icurly  at  the  greater  distance. 

In  conclusion,  we  feel  that  our  model  adequately  describes  the 
propagation  of  plane  sound  waves  in  a real  pipe  in  the  frequency  range 
of  at  least  0.5  to  3.5  kHz,  at  sound  pressure  levels  up  to  at  least 
160  dB.  The  application  of  the  model  to  outdoor  propagation  is  expected 
to  be  strcdght forward.  The  transformation  from  plane  to  spherical  or 

(9) 

cylindrical  waves  is  well  known.  Furthermore,  the  absence  of  wall 
effects  will  represent  a simplification. 
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